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Abstract: We present stellar yields calculated from detailed models of low and intermediate-mass 
asymptotic giant branch (AGB) stars. We evolve models with a range of mass from 1 to 6Mq, and 
initial metallicities from solar to 1 /200 th of the solar metallicity. Each model was evolved from the 
zero age main sequence to near the end of the thermally-pulsing AGB (TP- AGB) phase, and through 
all intermediate phases including the core He-flash for stars initially less massive than 2.5M0. For 
each mass and metallicity, we provide tables containing structural details of the stellar models during 
the TP-AGB phase, and tables of the stellar yields for 74 species from hydrogen through to sulphur, 
and for a small number of iron-group nuclei. All tables are available for download. Our results have 
many applications including use in population synthesis studies and the chemical evolution of galaxies 
and stellar systems, and for comparison to the composition of AGB and post- AGB stars and planetary 
nebulae. 

Keywords: stars: AGB and post-AGB — abundances — ISM: abundances — planetary nebulae: 
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1 Introduction 

Stars with initial masses in the range ~ 0.8 to 8M©, depending on the initial metallicity Z, will pass through 
the thermally-pulsing asymptotic giant branch (TP-AGB) phase before ending their lives as white dwarfs. The 
structure and evolution of low and intermediate mass sta rs prior to and during the AGB has been previously 
discussed in some detail fsee lBusso et al ] l!999l :lH erwidl2005l and references therein) . An AGB star is characterized 
by two nuclear burning shells, one burning helium (He) above a degenerate carbon-oxygen core and another burning 
hydrogen (H), below a deep convective envelope. In-between lies the intershell region composed mostly of 4 He. 
The helium-burning shell is thermally unstable and flashes or pulses every 10 4 years or so. After the occurrence of 
a thermal pulse (TP), mixing episodes may occur that bring the products of nuclear burning from deep inside the 
star to the stellar surface. These mixing events, or third dredge-up (TDU) episodes, bring the products of partial 
He-burning (mostly 4 He and 12 C) into the envelope. The TDU is the mechanism for turning (single) stars into 
carbon stars, where the C/O ratio exceeds unity in the surface layers. Following dredge-up the H-shell is re- ignited 
and the star enters a phase of quiescent H-burning known as the interpulse phase. The TP-AGB is defined as the 
phase from the first TP to the time when the star ejects its outer envelope, terminating the AGB phase. 

In the most massive AGB models (M > 4Mq, depending on the metallicity) hot bottom burning (HBB) can 
occur when the base of the convective envelope dips into the top of the H-shell resulting in a thin layer hot enough to 
sustain proton-capture nucleosynthesis. The convective turn-over time in the envelope is ~ 1 year that ensures that 
the whole envelope is exposed to the hot region a fe w thousand times per interpulse period. There is observational 
evide nce that HBB occurs in massive AGB stars ([Wood et all 119831 : ISmith fc Lambert] [l98 9: McSa venev et all 
120071 1 in the Large and Small Magellanic Clouds (LMC and SMC). HBB converts 12 C into 14 N and if the TDU is 
operating, HBB will prevent the atmosphere from becoming carbon rich (Boothroyd, Sackr nann. fc Ahernlll993h . 
The copious amounts of 1 4 N pr oduced in this case will be primarjQ owing to the primary 12 C being dredged from 
the He-shell. iFrost et al.l (1 19981 ) found intermediate-mass AGB evolution to be sensitive to the initial composition 
and the mass-loss law used in the calculation, with models becoming C-rich near the tip of the AGB if TDU 
continues after HBB has been shut off. The occurrenc e and duration of HBB depends upon the initial metallicity, 
mass loss and convection (I Ventura fc D'Anto na 2005a) where the minimum mass for HBB is pushed to lower mass 
in lower metallicity models ~l|Frost et al.l 199sh . 



Owing to fact that calculating a TP-AGB model is a computationally intensive task, synthetic AGB models, 
which use fitting formulae to model the evolution quickly, have proved to be a successful approach to population 
syntheses studies that require large stellar populations. Historically this approach was validated by the fact that 
the stellar luminosity on the AGB is nearly a linear function of the H-exhausted core mass (e.g. lPaczvnski|[l975l : 
iRenzini fc Volilll98ll : iGroeneweeen fc de Jondll993h although we know now that this relation breaks down for 

produced from the H and 4 He initially present in the star. 
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stars undergoing HBB (Bloec ker fc Schoenbe£ne^ |lj9iyLattarrzk |l992|). Synthetic AGB mo dels have successfully 
been u sed to model AGB populations dGroenewegen fc de Jong _ 19931; iMarigo et all 1 19961 ) and compute stellar 
yields (|van den Hoek fc Groenewegen|[l997l : iMarigd 200ll : llzzard et al.ll2004l ). Many of the parameterizations used 
in synthetic evolution studies are derived from detailed stellar models, such as the growth of the H-exhausted core 
with time, and as such are only a ccurate over the range in mass and me tallicity of the stellar models they are based 
upon. An example is provided bv lvan den Hoek fc Groenewe gen ( 1993) who compute AGB yields for initial masses 
between 0.9 and 8M0 whereas the lBoot4irovd"^ri5ack*mannl (119881 ) interpulse-period-core mass relation they use 
was only derived for stars with initial masses between 1 and 3M© . What affect this has on the yields is unclear since 
this relation will affect the number of TPs during the TP- AGB phase and hence the level of chemical enrichment. 
Dredge-up and HBB affect these relationships but owing to our lack of understanding about these phenomena, 
there are still considerable uncertainties that affect the detailed computations. Recent impro vements in computer 
power mean that grids of detailed AGB models can now be produced in a reasonable time ijVentura et al ] |2002l : 
iKarakas et al.ll2002l ; iHerw ig 2004b) although producing yields from N > 20 AGB stars for any given metallicity 
range is still challenging. For this reason synthetic models are still preferred for some applications. 

In an effort to try and quantify the contribution of AGB stars to the production of elements in the Galaxy 
we have calculated grids of detailed AGB models covering a substantial range in initial mass and metallicity. 
The aim of this paper is to provide the results of these grids in a format suitable for many applications. We 
first present the AGB stellar structure in a format useful for comparison to other detailed models, or for use in 
synthetic AGB algorithms. Second, we provide the stellar yields. These could be useful for different applications 
including chemical evolution studies, for comparison to the composition of planetary nebulae (PNe) and pre-solar 
grains. The data are provided in a tabulated format available for download via the internet, where the appropriate 
reference for these tables is this paper. 

These models were previously discussed in detail infjKarakas (2003) and in a number of other studies that help 
to hig hlight the importanc e of these models for both de t ailed and synthetic AGB calculati o ns. These studies in- 



clude IKarakas et all (l2002l)JlKarakas fc Lattanzid (120031) . Ilzzard et all (l2004l ). lLugaro et all l|2004l ). IKarakas et all 
(|2006bl ). lLugaro et all (|2007l ). Ilzzard et al.l (|2007T l and lKarakas et al.l l|2007f ). These later studies used the structure 
models presented here and changed the reaction rates and/or initial abundance s to study particula r problems e.g. 
the effect of reaction rate uncertainties on the production of 19 F in AGB stars (jLugaro et al.1120041 ). 
We begin with a brief description of the numerical method. 



2 The Numerical Method 

The n umerical method we use has been previously described in detail (e.g. iLugaro et al. 1 12004 IKarakas et alj 
l2006bl ). Here we summarize the main points relevant for this study. We compute the structure first and then 
perform detailed nucleosynthesis calculations. 

2.1 The Stellar Structure models 

The stellar st ructure models we r e calc ulated with the Monash version of the Mount Stromlo Stellar Structure 
Program; see iFrost fc Lattanzid (|1996| ) and referen ces therein for details. Mass loss on the first giant branch is 
included using the Reimer's mass-loss prescri ption (lReimerslll975 . here after R75) with the parameter r\ = 0.4; 
on the AGB we use the formulation given by IVassiliadis fc Wood! l| 19931 . hereafter VW93) in all models unless 
indicated otherwise. We calculate two models (5Mq, Z = 0.02 and Z = 0.0001) using R e imer' s mass loss on the 
AGB with the parameter r\ = 3.5. We used the solar abundances of lAnders fc Grevessd (Il989l) for the Z = 0.02 
models. For the Z = 0.008 and 0.004 models we used C, N, and O abundances from lRussell fc Dopital (|1992| ) that 
are appropriate for stars in the Large and Small Magellanic Clouds where many AGB studies are performed. We 
assumed scaled-solar abundances for the Z — 0.0001 models. 

All models were calculated from the zero-age main sequence to near the end of the TP-AGB phase. For 
models with initial masses below 2.5Mq we also evolve the models through the core He-flash. While this is an 
approximation to the true 3D nature of the core flash ((Dearborn et all 12006). it at least allows us to model the 
entire evolution of the star self-consistently wit hout resorting to using zero-age horizontal branch models or other 
techniques to get around this event (see for e.g. Stancliffc et al. 20q3). 

T he occurrence of the TDU and HBB depen d critically upon the convection model used l|Frost fc Lattanziol 

119961 ; lMowlavilll999l ; IVentura fc D'Antonall2005al ) and the method for determining convective borders. We use the 
standard mixing-length theory (MLT) for convective regions with a mixing-length parameter a = l/Hp — 1.75, 
and determine the border by applying the Schwarzschild criterion. Hence we do not include convective overshoot, 
in the usual sen se. Instead we search for a neutral border to the convective zone in the manner described by 
lLattanzid (|1986| ). This method has be en shown to increase the efficiency of the TDU compared to models that 
strictly use the Schwarzschild criterion l|Frost fc Lattanziolll996l ). 

We will not provide a detailed discussion of t he un c ertainties that affect our models but we refer the i nterest ed 
reader to lFrost fc Lattanzid (|l996l ). lBusso et all (|l999l ), lGorielv fc Mowlavil pOOOl ). I Ventura fc D'Antonal (|2005al lbh 
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Table 1: Grids of stellar masses for each Z, noting if the models experience the core He-flash (CHe), the 
third dredge- up (TDU), and hot bottom burning (HBB). 
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Z = U.UUo 


Z = U.UU4 


Z = 1U 


1.0 


CHe 


CHe 


CHe 
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TDU 
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TDU 


TDU 


TDU 
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TDU 


TDU 


TDU 


TDU 
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TDU 


TDU 


TDU 


TDU,HBB 
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TDU 


TDU 


TDU 


TDU,HBB 
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TDU 


TDU, HBB 


TDU, HBB 


TDU, HBB 
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TDU, HBB 


TDUHBB 


TDU, HBB 


TDU,HBB 
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TDU, HBB 


TDU, HBB 


TDU, HBB 


TDU,HBB 


6.5 


TDU, HBB 








7.0 








TDU,HBB CI 



a This model goes through central C-burning and is a super- AGB star. 



andl crwi g (2005) noting that this is a non-exhaustive list of publications dedicated to discussing uncertainties in 
the modelling of AGB stars. 

2.2 Post-processing Nucleosynthesis models 

We performed detailed nucleosynthesis calculations u sing a post-pr ocessing code that includes 74 species and 
time-dependen t diffusive mix i ng in all convective zones (|Cannonll 19931 ). The details of the nucleosynthesis network 
are outlined in iLuearo et all ([2004 ) but we remind the reader that we include 59 light nuclei and 14 iron-group 
species. We also add the fictional partic le g to count the number of neutron captures occurring beyond 62 Ni 
(jLattanzio et al.lll996l ; lLugaro et al.l(2004l '). 

The bulk of the 506 reaction rates are from the REACLIB data tables (Thielem ann et al1ll986h . based on the 
1991 updated version. Some of the proton, a and neut ron capture reaction rates have been updated according to 
the latest experimental results, see lLuearo et al.1 (|2004l ) for details. 

3 The Stellar Models 

The mass grid used for each metallicity is shown in Table [T] where we note if the model experienced the core 
He-flash, TDU and/or HBB. For each model AGB star we present structural information in tabulated form and 
these tables are available for download from the website http://www.mso.anu.edu.au/~akarakas/model_data/ 
Samples are given in Tables [5] and [3] Data is provided as a function of the TP number as a proxy for time. Each 
table contains the pulse number, and for that pulse we include the core mass M core , the maximum extent of the 
pulse-driven convection region M csri , the maximum duration of the convective pocket t csri , the amount of matter 
dredged into the envelope AMdredge, the parameters A and Xdu P (see below), and the maximum temperature in 
the He-shell ThosIicII (measured from the middle of the He-shell which is hotter than the bottom). We include the 
maximum temperature at the base of the convective envelope Tb ce , and the maximum temperature in the H-shell 
Thaheii; these were computed over the previous interpulse period. We then include the interpulse period, the total 
mass (measured at the current TP) M to t, the maximum radiated luminosity during the previous interpulse period 
Max L, and the maximum He-luminosity during the current TP Max LHe. All units are in solar units with the 
exception of temperatures, that are in kelvin, and all times which are in years. 

The TDU efficiency parameter, A, is usually defined according to A = AMdredge/AMh, where AMh is the 
amount by which the core mass has gro wn between the present an d previous TPs. In our tables we provide A 
along with the parameter Adup, defined bvlGorie lv fc M owlavil ()2000l ) to be AMdredge /M cs h. This quantity is useful 
for comparison to the results of lGorielv fc Mowlavil ((2000J) and also because it pr ovides a measure of the am ount 
of pulse-driven convective mass mixed into the outer envelope. In the study bv lGorielv fc Mowlavil (2000) they 
assumed a constant value of Ad up = 0.1 whereas we observe this parameter to increased as a function of TP number 
in models that experience the TDU. 
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We had originally aimed to evolve each model A log T e g ~ 0.3 off the AGB track, to the point where the entire 
outer envelope of the model had been lost. Owing to convergence difficulties near the tip of the AGB this was only 
achieved in one case: the 1M©, Z = 0.02 model where mass loss removed the entire outer envelope, leaving the 
model with a final remnant mass of O.573M0 and on the white-dwarf cooling track with a final logT c ff = 4.118. 
The low-mass models (M < 2M0) lost most of their outer envelopes, even if they did not leave the AGB track 
completely according to the definition above. For most of these models the remaining envelope mass was small, 
much less than O.IM0 and we would not expect further TPs on the AGB. For the more massive models the 
remaining envelope mass was large enough that further TPs would probably occur if we could evolve the models 
further along the AGB track. In all cases however the luminosity-driven superwind phase had started and the 
envelope mass was being reduced at a rapid rate, at a few x 10 Mq year -1 . HBB had ceased to operate in 
all models; this occurred once the envelope mass was reduced below about ~ 1.5M0. In Sj4]we discuss how we 
estimated the number of remaining TPs and the contribution of these pulses to the chemical enrichment of the 
model star. 

The convergence difficulties mostly occurred during the start of the final TDU phase, after the TP has peaked 
in He-luminosity. The reasons for the convergence problems are unclear but they are not associated with small 
values of the gas pressure to total pressure, j3, since previous dredge-up episodes occur without problem for similar 
values of /3. The convergence problems are sometimes associated with a rapid drop in the temperature at the base 
of the convective envelope and an increase in the stellar radius. Higher metallicity models, at a given initial mass, 
experience convergence difficulties earlier (i.e. at a larger envelope mass) than models of lower metallicity. This 
may indicate that the input physics, notably the opacities, are related to the problem. A thorough investigation 
into the cause of the convergence difficulties is required. 

Two models require a special mention and these are the 6 and 7Mq, Z — 0.0001 models. The 6M0 model was 
evolved through more than 100 TPs but by the time we ended the computation, little envelope mass had been 
removed, the superwind phase had not been reached and HBB was actively changing the com position of the outer 
envelo pe. For this model we estimate that another 100 could possibly occur when using the I Vassiliadis fc Wood! 
( 1993) mass-loss prescription. We do not at tempt to estimate the contribution from the remaining TPs owing to 
the complication caused by active HBB; see Izz ard et al.l {2006) for the difficulties associated with this. The 7M0 
model went through degenerate central carbon burning before ascending the TP-AGB phase. For this model we 
only computed 12 TPs before we ended the computation. We point out that it is unclear if this model would 
end its life as a massive ONe white dwarf or it would reach the Chandrasekhar limiting mass and explode as a 
low-mass Type II supernova. At the end of the computation, the final core mass was 1.126M0 and with inefficient 
TDU and a low mass-loss rate we suspect it would explode as a low-mass Type II supernova but further study is 
required to determine the final fate of this model. 

In Tables [5] and [3] we present the results for the 3M0, Z = 0.02 model, to show an example of the tables 
available. Note that these data are given as one table for each mass and Z in the on-line version. For the case 
shown the TDU begins at the 10 th TP when the core mass reaches O.62M0 . Adding up the AM^dge for each pulse 
gives a total of O.O787M0 core material that is added to the envelope over the course of the TP-AGB phase. This 
model does not experience HBB and that can be seen by inspection of Table [3] where the maximum temperature 
at the base of the convective envelope during the interpulse does not exceed 10 7 K. Other interesting comparisons 
can be made by considering the evolution of the interpulse period as a function of the core mass. From the data 
presented in Table [3] we observe an initial increase followed by slow d ecrease. We can compare to the core mass- 
interpulse period relations given in Boothroyd & Sackmann] (| 19881 ) or I Wagenhuber fc Groenewegenl (|l998f ) at the 
final core mass of the model, O.68M0. The first relation gives a value of 3.47 xlO 4 years and the second 3.73 x 10 4 
years. Both of these values are significantly smaller than the value that comes out of the detailed computation of 
5.47 x 10 4 years. Hence synthetic AGB simulations using either of these relations would result in many more TPs 
during the TP-AGB than we predict and hence more TDU episodes and a greater level of chemical enrichment. 



3.1 The Third Dredge- up and HBB 

Owing mostly to the operation of TDU episodes AGB stars are important for enriching the composition of the 
interstellar medium (ISM). Observat ional evidence suggests that the TDU occurs in stars as low as ~ 1.2M@ 
at solar metallicity and in the LMC ijFrogel et al.lll990lTWallerstein fc Knapdll998T ) and at ~ 1M in the SMC 
(jFrogel et alJll990f ). An examination of the tables presented above show that the minimum mass for the TDU 
at solar metalli city is 2.25M0 , and I.5M0 at the metallicities appropriate for the LMC (Z = 0.008) and SMC 
(Z = 0.004). In lKarakas et al] (|2002h we discussed in detail how these models do not predict eno ugh TDU at small 
enoug h core masses to account for the observations. The recent analysis of infra-red data by iGuandalini et ahl 
(2006) suggests that the minimum mass of stars experiencing the TDU has previ ously been greatly unde r-estimated, 
owing to incorrect estimates of bolometric magnitudes. However if the study of Stancli ffe et al.l l|2005h is anything 
to go on there is still a requirement for more efficient TDU in low-mass low-metallicity stars than we predict. 

Massive AGB stars that expe rience HBB are one of the favoured s ites of the abundance anomalies observed 
in galactic globular cluster stars l|Yong et al 1 l2003l ; lGratton et al.ll200l^ even if problems persist in comparisons 
betw e en the model predictions and the observed abund ance trends (|Denissenkov fc H erwig 2003; Fenner et al.l 
120041 : IVentura fc D'Antonall2005bl : iKarakas et al . 2006a). Observations of stars experiencing HBB have mostly 
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come from the LMC and SMC (IWood et al.lll983l : lPlez et al.lll993l ). with some observations of stars in our Galaxy 
( Gar cia-Hernandez et al ] |2006l . l2007t) ~ These observations indicate that most luminous AGB stars have an O-rich 
comp osition (where C/O < 1) with a low 12 C/ 13 C ratio, and that many are also rich in Li l|Smith fc Lambert! 
1990), and in elements produced by the slow-neutron capture process. Observations bv iMcSavenev et al l (|2007h 
have shown that two bright O-rich AGB stars are rich in nitrogen and depleted in carbon. These observations are 
the first observational confirmation of the long-predicted production of primary nitrogen by the combination of 
TDU and HBB in intermediate-mass AGB stars. 



4 The Stellar Yields 



The stellar yields can be downloaded from the website: http : //www.mso . aim. edu. au/~akarakas/stellar .yields/ 

We present the stellar yields in two ways. The first method presents the yields and abundances in the wind in- 
tegrated over the entire stellar lifetime and these are suitable for use in e.g. galactic chemical evolution studies. 
An example is given in Table [5] Each of the tables contain the following: 1) the nuclear species i, 2) the atomic 
mass A(i), 3) the nett stellar yield, defined below, 4) the amount of species i in the wind lost from the star 
mass(i)i ost (Mq), 5) the amount of i that would initially have been present in the wind mass(i)o (Mq). The 
quantity mass(i)o is simply the mass expelled during the stellar lifetime multiplied by the initial abundance. We 
next include 6) the average mass fraction of i in the wind (X(i)), 7) the initial mass fraction Xo(i) and 8) the 
production factor / denned by log[(A(i))/Ao(i)]. Note that in Table [5]the production factor has simply been 
labelled / whereas the definition is provided in the header of the on-line tables. Thi s last quantity / is useful for 
comparison to authors who present their results in a similar way (see Table 4 from iHerwi g 2004b) as well as for 
comparison to derived abundances from stars, where for e.g. an / ~ 0.2 would indicate an increase of that species 
by 0.2 dex at the surface. 

The definition of the yield that we use is given by the following expression 

M = £ [X(i)-X (k)]^-dt, (1) 

where Mj is the yield of species i (in solar masses), dM/dt is the current mass-loss rate, X(i) and Xo(i) refer to 
the current and initial mass fraction of species i, and r is the total lifetime of the stellar model. The yield can 
be negative, in the case where the element is destroyed, and positive if it is produced. We also present the total 
amount of i (in Mq) expelled into the ISM, noting that this value is always positive. 

It was previously discussed that some of the models do not lose their entire outer envelopes during the TP- 
AGB evolution owing to convergence difficulties near the tip of the AGB. This may be important for the stellar 
yields because, depending on the amount of remaining envelope, these models may experience further TPs and 
TDU episodes. Even one more TDU episode could significantly alter the surface abundances at this stage because 
the envelope mass is small leading to less dilution. The affect on the stellar yields is probably smaller because this 
quantity is integrated over the entire stellar lifetime. 

For each model we estimate the number of remaining TPs; this is shown in Table [3] along with the TDU 
efficiency parameter from the last computed TP, Af. Models not included in this table have small enough envelope 
masses that we estimate that no further TPs would take place. For the models in Table H] we present two sets 
of stellar yields, one with no contribution from these extra pulses; this is where we remove the envelope at its 
given composition and assume that A = during remaining pulses. The other set of yields were computed using 
a constant value of A for the remaining TPs. This latter set are our "standard" set and we recommend these be 
used. The set with A = are useful for exploring model uncertainties. 

From inspection of Table [4] an example of a model where remaining TPs might affect the yields is the 6M© , 
Z = 0.02 model where we estimate that 5 more TPs could occur. The magnitude of the change to the surface 
composition from these final TPs depends on one big unknown: The behaviour of A wit h decreasing envelope 
mass. Detailed stel lar models suggests that A decreas es with decreasing env elope mass (jStraniero et al ] ll997l : 
iKarakas et al.l |2002| . and data presented here) but IStancliffe fc Jeffervl l|2007h find results to the contrary and 
instead suggest that there is no reduction in A, at least for the 1.5Mq models discussed in that paper. 

In Table[5]we present a set of stellar yields for the 3Mq, Z — 0.02 model. From Table^]we estimate that one 
more TP would occur for this model and the yields presented here include this final TP. In Table [5] we only show 
results for the 39 nuclear species that had stellar yields with mean magnitudes greater than ~ 10 -10 and hence 
neglect many short-lived unstable isotopes that are included in the nucleosynthesis network as well as in the tables 
that are available for download. Using this criterion we have included results for the important radio-nuclides 26 Al 
and 60 Fe, noting that we set the production factors to be zero for these models owing to an initial zero abundance. 

For the 3Mq, Z = 0.02 model we estimate that the last TP will occur when the total stellar mass was ~ 1.35M© 
and the envelope ~ 0.66Mq. From inspection of Table [3] it is clear that this last TDU episode occurs at a much 
smaller envelope mass than all previous mixing episodes, where M onv > 2Mq. Consequently the change to the 
surface abundances was larger and this can be seen in Fig. [T] where we show the evolution of the C/O ratio at 
the surface. To see what effect this last TDU episode has on the stellar yields we compare the nett yield of 12 C. 
Without this final TDU episode the nett yield of 12 C is 1.354x10 _2 Mq whereas with that final mixing event the 
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Figure 1: The evolution of the C/O ratio at the surface during the TP-AGB phase of the 3M Q , Z = 0.02. 
The final C/O ratio for this model is 1.43. 



yield is 1.496x 1O" 2 M , an increase of about 11%. On the other hand the change to the final surface abundance 
of 12 C is about 20%. 

Before we move on it is important to discuss the effect of the initial mass on the degree of chemical enrichment. 
It was mentioned above that the 6M0 Z = 0.02 model had an extra 5 TPs left, given the final envelope mass of 
1.53M0 from the end of the detailed calculation. While it would seem obvious that a greater number of remaining 
TPs would lead to a larger enrichment of the envelope, we need to remember that the He-intershell is about a 
factor of 10 smaller in an intermediate-mass AGB star. For this reason the amount of matter dredged into the 
envelope at each TDU episode is smaller, regardless of the value of A. This can be seem from inspection of the 
tables given in Sj3] where the maximum mass of the pulse-driven pocket (approximately equal to the mass of the 
intershell at the TP) is O.OI6M0 for the 3M0, Z — 0.02 model whereas this value is reduced to O.OO33M0 in the 
6M0, Z = 0.02 case. The mass of the envelope in the 6M0 model was also larger at each TDU episode resulting in 
more dilution of dredged-up material. For this massive AGB model, the final surface abundance of 12 C increased 
by about 26% as a consequence of the 5 extra TPs, only marginally more than the 20% increase from only one 
TP for the 3M Q , Z = 0.02 CQ/SG, ELS discussed above. 

The stellar yields that we present here may be slightly different to those presented in iKara kas (2003]) because 
we have improved the synthetic AGB algorithm used to estimate the surface enrichment from the final TPs. The 
main improvement is taking intershell compositions for the last TP for all 74 species from the detailed stellar 
models. Previously we estimated the intershell compositions for about 15 species (e.g. 4 He, 12 C, ie O, 22 Ne etc) 
and assumed the rest did not change. We also did not have intershell compositions for each mass and Z, so in 
some cases we assumed that compositions appropriate for e.g. the 3M Q , Z = 0.004 model were appropriate for 
the 2.5 and 3.5M0 models of the same Z. We have also improved the input data required for the synthetic code, 
e.g. the final computed interpulse period and TDU efficiency. These improvements change the final remnant mass 
(at least slightly) and this also results in slight changes to the stellar yields. 

4.1 Planetary nebulae 

After the outer envelope has been ejected t he stars evolve thr ough the brief post- AGB and planetary nebula (PN) 
phases of the Hertzsprung- Russell diagram l|van Winckei2003T ) . The glowing gas of the nebula is the remnant of the 
outer convective envelope of the progenitor AGB star and hence its composition should reflect the nucleosynthesis 
that took place during the final few TPs on the TP-AGB. 

Planetary nebulae are important because their compositions can be determined fairly accurately (at least for 
some species); because they are bright, abundances can be derived from PNe as far away as M31, and because the 
gas is ionized by UV photons we obtain information about elements normally very difficult to study, in particular 
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the noble eases He, Ne, Ar, as well as C, N, O ijKaledl 19781 : iHenrvl 19891 : iDopita et al.ll 19971 : [Stanghellini et all2 000). 
Abundances of heavy elements synthesized by the s p rocess including Ge, Se, Kr, Xe Br and Ba l|Sterling et al.l 
120021 : ISterling fc Dinersteinll2007( : ISharpee et ai1l2007n are also observed to be enhanced in some nebulae. Owing 
to the importance of PNe abundances we have included stellar yields for PNe in particular. These yields are the 
abundance by number in the matter lost from the star from the final two TPs and as such are good for comparison 
to the composition of the nebular gas for PNe in the Galaxy, LMC and SMC. We have included PN yields for the 
low-metallicity Z — 0.0001 models but few (if any) PN around today are expected to have evolved from such a 
low-metallicity population of stars. PNe in the Milk y Way Galaxy are predom inantly found in the disk with oxygen 
abundancefl in the range 7.73 < loge(O) < 9.09 (jStanghellini et all 120061 ) compared to the solar value of 8.66 
(|Asplund et al.l l2005). Only a few PN have been found in the halo and there are only a handful (JV ^ 12) known 
that are defined to be metal-poor, with oxygen abundances less than a factor of 4 below solar (jDinerstein et all 
l2003ft . 

In Table [H] we show an example of the yields available for a PN originating from a 3M© progenitor of solar 
composition. As for the yields presented previously we compute two sets for models with synthetic TPs and for 
this case we show the PN set with the final synthetic TP included. The first entry in each on-line table is the 
amount of matter expelled during the last two TPs, here we include this value in the table caption. The columns 
contain the species i, the atomic weight A(i), followed by the nett yield, the mass lost, the initial mass fraction, 
and the average mass fraction integrated over the final two TPs. The next quantity is the abundance (by number) 
of the matter lost in the wind over the final two TPs, compared to the number of hydrogen atoms, Ni/Nn- This 
quantity is often used in the literature in the form log e = log(Ni/Nn) + 12. 

Comparing results for 12 C from Table[S]to Table[5]we see that the nett yield is lower in the PN table by ~ 27%, 
this is consequence of less matter being lost in the final two pulses than over the whole lifetime, even though the 
average 12 C abundance is higher by ~ 8%. Remember that our definition of a stellar yield defined in Equation [T] 
essentially subtracts the amount of i in the wind lost from the star (i.e. » M^t* 1 x by the amount of i 

that would have been initially present in that wind (i.e. M^ 3,1 x Xo(i)). Note that the 3M Q , Z = 0.02 model lost 
1.5O8M0 over the last two TPs compared to 3.0 — M{ = 2.318M© over the entire lifetime (see TableH]for the final 
mass) . 

That the yield difference was not larger is a consequence of the VW93 mass-loss prescription. This prescription 
ensures that the mass-loss rate is low until the superwind begins at which stage most of the envelope is lost in a 
few TPs. For the 5M©, Z — 0.02 model we computed an evolutionary sequence with the R75 mass-loss law with 
r] — 3.5. For this model the nett PN yield of 12 C is 3.42 times lower than the total yield, compared to 1.68 times 
lower for the VW93 case. The total amount of matter expelled is O.75M0 and O.58M0, for the VW93 and R75 
cases, respectively. Note that the average 12 C abundance in the wind expelled over the last two TPs increased in 
both models, where the R75 model saw a 66% increase compared to 63% for the VW93 case. 

4.2 Limitations of the current models 

We discussed in the previous section that we did not include a partial mixing zone at the deepest extent of each 
TDU episode to produce a 13 C pocket. This can be considered a limitation since the neutrons produced by the 
13 C(a, n) ls O reaction can affect the composition of the He-shell. However the light elements (hydrogen through 
to phosphorous) are generally not significantly altered by neutron captures, at least at solar metallicity. There 
are exceptions to this, including th e nucleosynthesis of 19 F that requires th e operation of the 13 C neutron source 
in low-mass AGB stars of ~ 3M Q y orissen et all 1 19921 : iLugaro et al]|2004D . Sodium, the magnesium and silicon 
isotopes, and 31 P are affected by neutron captures during the TP where neutrons are released by the 22 Ne source 
at higher neutron densities (~ 10 10 cm -3 ), noting we include this in our computations. Neutron captures will 
impact more significan tly on light-element nucleosynthesis in more massive or lower metallicity models, as pointed 
out bv lHerwidl2004al lbl). 

One limitation of our models, at least in terms of their use for galactic chemical evolution studies, are the 
initial C, N and O abundances used for the Z — 0.008 and 0.004 models. These are not scaled solar but reflect 
the composition of the LMC and SMC, respectively. Our models are therefore very useful for comparing to the 
composition of AGB stars and PN in the Large and Small Magellanic Clouds, which was their intended use, but 
one must be careful when using the C, N and O yields for comparison to populations that have a similar [Fe/H] 
content of the Clouds but a different chemical history. We have started to compute yields assuming a scaled-solar 
composition for those metallicities and these will be available in a future publication. 

Another limitation is that our lowest metallicity models are Z = 10~ 4 whereas there are many new discoveries 
of low-mass halo stars with [Fe/H] < -3 (|Beers fc Christliebll2005h , some of which could have been polluted by 
a former AGB companion. We do not intend to co mpute AGB mod els with such low metallicities in the future 
but we refer the reader to Simon Campbell's models ijCamp bell 2007). Besides that 7M , Z = 0.0001 super- AGB 
star we do not compute models for any more of these objects. Yields for these stars should be available in the 
future (Doherty, Lattanzio & Siess, in prep aration), no ting that models of these stars are proving to be even more 
challenging than for CO-core AGB models (Sicss 2006). 



2 used as a proxy for metallicities because Fe is difficult to measure owing to condensation into dust. 
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We do not include the effect of extra-mixing processes into our computations, and by this we mean mixing 
events beyond those predicted to occur by standard models such as ours e.g. the first and second dredge-up, the 
TDU and HBB. Observations suggest that extra-mixing pr ocesses reduce the 12 C/ 13 C ratio in low-mass (m < 2Mq) 
red giants stars below that predicted by standard models (jCharbonnell 19941 ) . Observations of A GB stars and pre- 
molar grain analysis indic a te that a similar type of extra-mixing occurs in low-mass AGB stars (jBoothrovd et all 
Il995l : lAbia fc Iser4|l997l : iBusso et al"1ll999l : iNollett et all 120031 ) While the physica l mechanism for the extra- 
mixing is unknowr0 it is commonl y thought to be rotat ion (e.g. iPalacios et"a ,l. 2006), alth ough we a l so ref er to 
the mixing recently discov ered by (iEggleton et al.l l2006h and investigated further by lEggleton et al.l l|2007f ) and 
(jCharbonnel fc Zahdl2007h . 

4.3 Comparison with other yields 

In thi s section we provide a brief comparison of our yields to other authors. iKa rakas ( 200^) and llzzard et al.l 
(I2004T ) provided a detailed comparison of Izzard's synthetic yields to our detaile d models (which his synthetic 
AGB algorithm is based upon) , and to those from other synthetic AGB models ijvan den Hoek fc Groenewegenl 
Il997l : iMarigo 200l]). We do not repeat that discussion here but instead we present diagrams in Figs. l2lto llll and 
we discuss results for 12 C and 14 N. In Fig. |3]we show the stellar yields of 12 C, a nuclei representative of low-mass 
AGB evolution where the effects of TDU episodes are important. In Fig. [6] we show the yields of 14 N, a nuclei 
that tracks HBB nucleosynthesis. In each figure we show the results for three metallicities, Z — 0.02, 0.008 and 
0.004 compared to yields from synthetic AGB models, se e the figure captions for more details. In each figure we 
have weighted the yields by the three-component IMF of iKroupa et al.l (|l993l ) such that 

Ift = |^ = «Mo)Mi ) (2) 

where j/j is the weighted yield and Mi is given by equation [T] 

From the Figs. [4] and |B] we see that our yields are almost the same as llzzard et alJ (|2004l ) owing to the f act that 
these synthetic models were tweaked to match our results. Our yields are similar in behaviour to those of iMarigol 
(2001), although because her models experience deeper TDU at a lower core mass, the yields of He-sh ell material 
e.g. 4 He, 12 C are higher. We notice significant differences with Ivan den Hoek fc Groenewegenl (jl997t l especially 
m reg ards to 14 N. This is owing to their simplistic treatment of HBB nucleosynthes i s whi ch under-predicts the 
amount of CNO cycling compared to the other computations. IStancliffe fc Jeffervl |2007) provided a detailed 
comparison between yields from 1.5Mq models computed with the STARS code to our yields, with the summary 
that their models experience deeper TDU and hence go in a similar direction to Marigo's, that is larger yields of 
material that is dredged from the He-burning shell. 

In Fig. [5] we show the (unweighted) yields of 15 N and 19 F as a function of the initial mass and metallicity. In 
this figure we include the Z — 10~ 4 to show the effect of lower metallicity stellar evolution on the yields. The 
maximum 19 F yield is pushed to lower mass, owing to a higher He-burning shell and deeper TDU. The yields of 
15 N are negative, indicating this isotope is destroyed, and scale with metallicity except in the intermediate-mass 
Z = 10~ 4 models, where we find that 15 N is produced. This is as a result of breakout from the 14 N(p, 7) 15 
reaction at the temperatures and densities of HBB in these lower metallicity models. The production of 15 N at 
lower metallicities may have important consequences for the evolution of 15 N in the interstellar medium, especially 
since AGB stars are not considered producers of this rare isotope . 

We compare the Z = 10" 4 models to the AGB models of iHerwid (|2004bh of the same metallicity. This 
comparison is made easier by the fact that both sets of models assumed a sc aled-sola r comp osition. Different 
mass loss was assumed where we used VW93 mass loss and Herwig applied the iBlockerl (11995T ) prescription with 
the parameter r/B = 0.1, that resulted in more rapid mass loss than we observed. This can be demonstrated by 
considering the 2Mq, Z — 10 -4 model as an example. Our 2M© calculation experienced 26 TPs compared to 
7 for Herwig's model. Even though the maximum amount of matter dredged into the envelope at any one TP 
was almost the same (1.112xl0~ 2 compared with Herwig's 1.171X 10 _2 Mq), the total amount of material dredged 
into the envelope over the whole TP- AGB was much larger in our case (0.2183Mq) compared to Herwig's model 
(5.038xlO~ 2 M Q ). For this reason our yields and average surface mass fractions are, in general, larger. The average 
abundance of 12 C in the wind of our model was ~ 2 times larger, of 22 Ne 13 times larger and of 23 Na 40 times 
larger. However the average abundance of le O was about 8 times larger in Herwig's model. These differences can 
be attributed to 1) more TPs in our model, leading to more TDU episodes and hotter nucleosynthesis conditions 
in the He-shell, and 2) different intershell compositions caused by diffusive convective overshoot used in the Herwig 
model. Overshoot applied to the border at the base of the convective pocket during a TP mixes some CO core 
matter into the intershell, resulting in increased 12 C and le O abundances. These yields differences are similar for 
all masses in the range 2 - 6Mq , where the magnitude of the difference increases with increasing mass owing to 
the greater number of TPs experienced by our intermediate mass models (4, 5 and 6Mq). 



3 and this is the main reason why we do not include it! Another reason being that it is essential to have standard models 
to compare with before testing new physics. 
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Figure 2: Weighted yield of 1 H as a function of the initial mass for the Z — 0.02 (top), the Z = 0.008 
(middle) and the Z = 0.004 models ( bottom). We show result s from our calculations (black s olid points), 



Ivan den Hoek fc Groenewegenl (|l997l ) (open magenta squares) jForestini fc Charbonnell (119971) (solid green 
square s), lMarigd "( 200 ll ) (open red circles). TVentura et alJ ( 20021 ) (solid aqua triangles) and llzzard et al.l 
(120041) (blue crosses). Forestini & Charbonnel (1997) and Ventura et al. (2002) do not provide yields for 
Z = 0.008 and cover a narrower mass range, between 3 and 6M Q . 
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Figure 3: Weighted yield of 4 He as a function of the initial mass for the Z = 0.02 (top), the Z — 0.008 
(middle) and the Z — 0.004 models (bottom). Symbols are the same as in Fig. [21 
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Figure 4: Weighted yield of 12 C as a function of the initial mass for the Z = 0.02 (top), the Z = 0.008 
(middle) and the Z = 0.004 models (bottom). Symbols are the same as in Fig. [21 
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Initial mass (M Q ) 

Figure 5: Weighted yield of 13 C as a function of the initial mass for the Z = 0.02 (top), the Z = 0.008 
(middle) and the Z — 0.004 models (bottom). Symbols are the same as in Fig.O 
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Figure 6: Weighted yield of 14 N as a function of the initial mass for the Z — 0.02 (top), the Z — 0.008 
(middle) and the Z — 0.004 models (bottom). Symbols are the same as in Fig. [21 
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Figure 7: Weighted yield of 16 as a function of the initial mass for the Z = 0.02 (top), the Z = 0.008 
(middle) and the Z — 0.004 models (bottom). Symbols are the same as in Fig. [31 
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Figure 8: Unweighted yield of 15 N (top panel) and 19 F (lower panel) as a function of the initial mass and 
metallicity Z. 
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Figure 9: Weighted 22 Ne (top panel) and 23 Na (lower panel) yields and for the Z = 0.02, the Z = 0.008 
and Z ~ 0.004 models. The yields of iForestini fc Charbonnel ( 19971 ) are also shown. 



www.pu blish . csiro. au /journals /pasa 



17 



CO 
I 

o 

X 
Lf) 



CO 

cn o - 
m 

CM 

H — 1^ 

O I 

o 

» in 
>- 



o - 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — T — I 1 1 1 1 1 1 1 1 [ 

. Z=0.02 / -----O 

□ Z=0.008 

O Z=0.004 / 

■ FC97 Z=0.02 O / _ 

a FC97 Z=0.005 'V /' jg - 




3 4 5 

Initial mass (M Q ) 



I 

O 



©r 



X 



CM 



CD 

y 



co 
I 

O 



I 

O 

X 



O - 



t — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — t — i — i — i — i — | — i — i — i — r- 

. Z=0.02 ® 
□ Z=0.008 / \ 

O Z=0.004 / • 

■ FC97 Z=0.02 • \ 

a FC97 Z=0.005 I \ 




Initial mass (N/U) 



Figure 10: Weighted 25 Mg (top panel ) and 26 Mg (lower panel) yield s and for the Z = 0.02, the Z = 0.008 
and Z — 0.004 models. The yields of iForestini fc CharbonneT ( 19971 ) are also shown. 
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Figure 11: Weighted 26 Al (top panel ) and 27 Al (lower panel) yields and for the Z = 0.02, the Z = 0.008 
and Z = 0.004 models. The yields of iForestini fc Charbonnel (|l997f ) are also shown. 
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5 Conclusions 

We present stellar structure data and yields from detailed AGB models with initial masses between 1 to 6M0, 
and with metallicities in the range from solar to l/200 th of the solar metallicity. All models were evolved from 
the zero-age main sequence to near the tip of the AGB track, from there we used a synthetic AGB algorithm to 
estimate the contribution of any remaining TPs. All results are presented in tabulated form and are available for 
download. These data are useful for many applications including comparison to other detailed AGB models, for 
use in galactic chemical evolution studies and for comparison to the composition of AGB stars, planetary nebulae 
and other post- AGB objects, as well as pre-solar grains. Such a comparison may help shed light on the model 
uncertainties that affect the yields, such as the mass-loss rate or nuclear reaction rates. Future work will involve 
running scaled-solar models for the intermediate-metallicities (Z = 0.008, 0.004) and adding more nuclei in the 
post-processing nuclear network to enable us to provide yields for s-process elements. Future work should also 
look at some of the reasons for convergence problems in the stellar models near the tip of the AGB, and include 
more physics into the calculations such as extra-mixing in low-mass AGB stars. 
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Table 2: Data for the 3M , Z = 0.02 model. 



Pulse 




-^^core 






M csh 






^csh 


AAfdrodgo 


A 


Adup 


1 


5. 


763726E- 


■01 


2 


237684E- 


-02 


3 


339341E+02 


000000E+00 


000000E+00 


000000E+00 

V'.V/V/V/V/V/V/l J | V / V / 


2 


5. 


782992E- 


-01 


2 


551377E- 


■02 


2 


.015471E+02 


000000E+00 


000000E+00 


000000E+00 

V/.V/V/V/V/V/V/l J | V / V / 


3 


5, 


.811490E- 


■01 


2 


586406E- 


-02 


2 


.061110E+02 


000000E+00 


000000E+00 


000000E+00 

V/.V/V/V/V/V/V/l J | V / V / 


4 


5. 


850938E- 


■01 


2 


.622330E- 


-02 


1 


793678E+02 


000000E+00 


000000E+00 


000000E+00 

V'.V/V/V/V/V/V/l J | V / V / 


5 


5, 


898175E- 


-01 


2 


.566504E- 


-02 


1 


.954225E+02 


0.000000E+00 


0.000000E+00 


0.000000E+00 


6 


5. 


952812E- 


-01 


2 


.515078E- 


-02 


1 


.755154E+02 


0.000000E+00 


0.000000E+00 


0.000000E+00 


7 


6. 


.012007E- 


-01 


2 


.422720E- 


-02 


2 


.020128E+02 


0.000000E+00 


0.000000E+00 


0.000000E+00 


8 


6. 


074721E- 


-01 


2 


.330387E- 


-02 


1 


.798662E+02 


0.000000E+00 


0.000000E+00 


0.000000E+00 


9 


6. 


139326E- 


-01 


2 


.226830E- 


-02 


1 


.702774E+02 


0.000000E+00 


0.000000E+00 


0.000000E+00 


10 


6. 


205002E- 


-01 


2 


.124983E- 


-02 


1 


.376081E+02 


2.955794E-04 


4.524140E-02 


1.390973E-02 


11 


6. 


269240E- 


-01 


2 


.039099E- 


-02 


1 


.702900E+02 


6.338358E-04 


9.432908E-02 


3.108412E-02 


12 


6. 


.331832E- 


-01 


1 


.966518E- 


-02 


1 


.570454E+02 


1.435995E-03 


2.083254E-01 


7.302221E-02 


13 


6. 


.389904E- 


-01 


1 


.907504E- 


-02 


1 


.575242E+02 


2.473533E-03 


3.414993E-01 


1.296738E-01 


14 


6. 


442509E- 


-01 


1 


.871735E- 


-02 


1 


.404630E+02 


3.429413E-03 


4.434169E-01 


1.832211E-01 


15 


6. 


.48971 1E- 


-01 


1 


.843846E- 


-02 


1 


.480161E+02 


4.516482E-03 


5.541952E-01 


2.449490E-01 


16 


6. 


.531649E- 


-01 


1 


.828820E- 


-02 


1 


.291233E+02 


4.803419E-03 


5.514702E-01 


2.626514E-01 


17 


6. 


571379E- 


-01 


1 


.772785E- 


-02 


1 


.474301E+02 


5.583823E-03 


6.362296E-01 


3.149746E-01 


18 


6. 


.607848E- 


-01 


1 


.762992E- 


-02 


1 


.466641E+02 


6.045818E-03 


6.549662E-01 


3.429294E-01 


19 


6. 


.641064E- 


-01 


1 


.734108E- 


-02 


1 


.509221E+02 


6.349206E-03 


6.777933E-01 


3.661368E-01 


20 


6. 


672819E- 


-01 


1 


.703131E- 


-02 


1 


.502906E+02 


6.745815E-03 


7.082441E-01 


3.960831E-01 


21 


6. 


.702475E- 


-01 


1 


.687807E- 


-02 


1 


.387330E+02 


7.370591E-03 


7.589639E-01 


4.366964E-01 


22 


6. 


.729269E- 


-01 


1 


.684576E- 


-02 


1 


.398080E+02 


7.480979E-03 


7.443761E-01 


4.440867E-01 


23 


6. 


.753748E- 


-01 


1 


.643026E- 


-02 


1 


.353630E+02 


7.710814E-03 


7.766045E-01 


4.693057E-01 


24 


6. 


777956E- 


-01 


1 


.632255E- 


-02 


1 


.362144E+02 


7.740200E-03 


7.639663E-01 


4.742027E-01 


25 


6. 


.800572E- 


-01 


1 


.598930E- 


-02 


1 


.370616E+02 


6.662667E-03 


6.661482E-01 


4.166952E-01 



Table 3: Data for the 3M fl , Z = 0.02 model continued. 



Pulse 




T 1 

- 1 Hcshcll 




71 

- 1 bcc 




she ll 


i ntprm 1 1 qa 




A/A. t 




MaxL 




Ma-xT/Hp 

IV LcLJL J_J J. 


1 


1 


O A A 1 1 ATT' i A O 

.840110E+08 


2 


.505874E+06 


4 


OOOP-1 A TT 1 i AT 

.832516E+07 


A AAAAAAT? i AA 

0.000000E+00 


2 


A O A A / 'AT7 1 i AA 

.989960E+00 


3 


OA^AAC T7 1 i AO 

.297096E+03 


2 


OA/ 1 /! AO TT 1 1 A /I 

.896448E+04 


2 


1 


A/ 1 T7 1 i AO 

.969207E+08 


2 


.684860E+06 


5 


.026387E+07 


PT A A 1 O A 1 I i A A 

5.221891E+04 


2 


AOAATATTi i AA 

.989959E+00 


3 


AOTr'O A T7 1 i AO 

.937624E+03 


1 


1 O O A A ATT i A r 

.183492E+05 


3 


2 


.086133E+08 


2 


.863514E+06 


5 


C\ (\ A O ATTTI i AT 

.204397E+07 


A TO A O T f 1 1 i A A 

6.784875E+04 


2 


AOAAP"OT71 i AA 

.989958E+00 


4 


ron ,1 T 1 , 1 i AO 

.585147E+03 


2 


AOTAP'/^TT i A P" 

.987256E+05 


4 


2 


1 nco /loir' i no 

.195342E+08 


3 


A01 /I O A T7 1 i Ai 1 


rr 

5 


O T 1 AOOT7 1 i AT 

.371028E+07 


O A01 1 r^L 1 i A/1 

8.0311ooE+04 


2 


AOAACCTT 1 i AA 

.989955E+00 


5 


A/ICTTAOTT 1 i AO 

.245723E+03 


7 


CACOOOT7' i APT 

.595382E+05 


5 


2 


.269594E+08 


3 


1 aa aaa i ri i aa 

.190444E+06 


5 


) (\ A A T A T7 1 i AT 

.494672E+07 


o a r a r att i a /i 

8.659560E+04 


2 


AOAAr"1 TTi i AA 

.989951E+00 


5 


TTA1 TTTTi i AO 

.779177E+03 


1 


01 AA/TT? i A/"» 

.319066E+06 


6 


2 


.342317E+08 


3 


.365202E+06 


5 


.611659E+07 


A f)Tf A"n i A /I 

8.875659E+04 


2 


AAAA A f T7S i AA 

.989945E+00 


6 


.321060E+03 


2 


A A 1 AAA"n i A 

.281208E+06 


7 


2 


.394346E+08 


3 


.526668E+06 


5 


.707H9E+07 


A irffn A f f in i a A 

8.778255E+04 


2 


AAAAO m I i A A 

.989937E+00 


6 


/7A/7W7/7F1 j r\ a 

.797777E+03 


3 


AAAA-l fTTl i A /~> 

.222215E+06 


8 


2 


.440928E+08 


3 


.692101E+06 


5 


l 7n/ 1 OAOT7' i AT 

.796298E+07 


o r 1 O A r O T7 i A /I 

8.513058E+04 


2 


A O A A A 1 T7 1 i AA 

.989924E+00 


7 


AAOTAP'TT 1 i AO 

.263765E+03 


4 


OAOAOATT i A/"» 

.398030E+06 


9 


2 


/I701 /"* O I 1 i A O 

.478168E+08 


3 


OAAAAOT7 1 i AA 

.866263E+06 


5 


OTr'OAATT i AT 

.875260E+07 


O 1 O /I O O A T7* i A /I 

8.134384E+04 


2 


AOAAAfTTi i AA 

.989906E+00 


7 


TA/^TTOTTi i AO 

.706778E+03 


5 


.538496E+06 


10 


2 


.513340E+08 


4 


A/IOOATTT 1 i A/? 

.048897E+06 


rr 

5 


A r A /I O O T7 1 i AT 

.950483E+07 


T TIOAOOTT 1 i A A 

7.718683E+04 


2 


AOAOTAT7 1 i AA 

.989879E+00 


8 


A A A O A ATT 1 i AO 





77/ 1 OAAT7 l i A/ 1 

.77o800E+0o 


n 


2 


p* p* 1 oro i ^ i a o 

.551253E+08 


4 


ni r ADATTi i aa 

.215980E+06 


6 


AOT'OAO T7 1 i AT 

.026202E+07 


T O A A A O TTT i A A 

7.320287E+04 


2 


A O A O 1 A T7 1 i AA 

.989840E+00 


8 


C O /I T1 /I T7 1 i A O 

.584714E+03 


8 


/I^^OA/IT? 1 A/^ 

.464204E+06 


12 


2 


.589790E+08 


4 


.410556E+06 


6 


.096342E+07 


6.977448E+04 


2 


.989780E+00 


9 


.015570E+03 


1 


.072477E+07 


13 


2 


.630434E+08 


4 


.610824E+06 


6 


.166960E+07 


6.747551E+04 


2 


.989682E+00 


9 


.461654E+03 


1 


.447887E+07 


14 


2 


.682954E+08 


4 


.819258E+06 


6 


.234783E+07 


6.636544E+04 


2 


.989504E+00 


9 


.909229E+03 


1 


.847517E+07 


15 


2 


.715982E+08 


5 


.026276E+06 


6 


.293235E+07 


6.538550E+04 


2 


.989177E+00 


1 


.032697E+04 


3 


.002030E+07 


16 


2 


.781746E+08 


5 


.236472E+06 


6 


.348534E+07 


6.524187E+04 


2 


.988530E+00 


1 


.073885E+04 


4 


.387040E+07 


17 


2 


.808189E+08 


5 


.397694E+06 


6 


.389052E+07 


6.336352E+04 


2 


.987350E+00 


1 


.108444E+04 


4 


.838104E+07 


18 


2. 


.832818E+08 


5 


.603128E+06 


6 


.432855E+07 


6.284853E+04 


2 


.985032E+00 


1 


.144775E+04 


5 


.528658E+07 


19 


2 


.867869E+08 


5 


.794352E+06 


6 


.465809E+07 


6.148713E+04 


2 


.980963E+00 


1 


.176355E+04 


8 


.176797E+07 


20 


2 


.898100E+08 


6 


.002216E+06 


6 


.496616E+07 


6.025162E+04 


2 


.973891E+00 


1 


.207047E+04 


9 


.513553E+07 


21 


2 


.926001E+08 


6 


.200394E+06 


6 


.525625E+07 


5.917455E+04 


2 


.962824E+00 


1 


.236732E+04 


1 


.163203E+08 


22 


2 


.953747E+08 


6 


.541253E+06 


6 


.554327E+07 


5.878434E+04 


2 


.951519E+00 


1 


.266948E+04 


1 


.445130E+08 


23 


2. 


.972316E+08 


6 


.751874E+06 


6 


.573163E+07 


5.686887E+04 


2 


.936666E+00 


1 


.291685E+04 


1 


.500224E+08 


24 


2. 


.999315E+08 


6 


.557310E+06 


6 


.592716E+07 


5.622341E+04 


2 


.845449E+00 


1 


.314644E+04 


1 


.755376E+08 


25 


3 


.021500E+08 


6 


.128534E+06 


6 


.598145E+07 


5.471787E+04 


2 


.211516E+00 


1 


.323496E+04 


1 


.658657E+08 
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Table 4: Data used to estimate the number of remaining TPs. We include the final total and core mass 
after the final TP, as determined by the synthetic AGB program, the A used for the final TDU episodes, 
the total number of TPs (including the synthetic TPs) an d the number of synthetic TPs. The VW93 
indicates the 5M Q , Z = 0.02 model that was computed with IVassiliadis fc Wood! ( 19931 ) mass loss on the 
AGB whereas the R75 indicates data for the model computed with Reimer's mass loss with rj = 3.5. 







Zj — 


-f\ no 
u. uz 






Mo 


ilif 


1\ I 

M cJ 


Af 


1NO. IrS 


oyn. irs 


2.5 


1.1388 


(\ O O f 1 (~\ 

0.6630 


0.4584 


26 


1 


3.0 


1.3580 


0.6819 


0.6661 


26 


1 


3.5 


1.2675 


0.7181 


0.8583 


22 


1 


4.0 


1.1151 


0.7916 


0.9306 


19 


2 


5.0 (VW93) 


1.2545 


0.8757 


0.9098 


27 


4 


O.U {txt o ) 


1 A/l ^8 


U.o i oO 


n one 


1 9 


r. 
■J 


O.U 


1 1 QSzL 


n Q9Q4 


n Q/L71 


40 


r. 



o.o 


u.yyio 


U.yozy 


n q con 


40 


/ 






Zj = 


U.UUo 






2.5 


0.9978 


0.6626 


0.7760 


28 


1 


3.0 


1 AA77 

1.0077 


0.6934 


0.8612 


on 


1 


O.O 


1/11 oo 


U. 1 OOz 


U.yOoU 


01 


1 

1 


4.0 


0.9716 


0.8367 


0.9730 


24 


2 


5.0 


1.1631 


0.8860 


0.9540 


61 


4 


6.0 


1.1408 


0.9480 


0.9195 


73 


5 






Z= 


0.004 






2.5 


1.3578 


0.6707 


0.8165 


29 


1 


3.5 


1.1274 


0.8138 


0.9720 


24 


1 


4.0 


1.0533 


0.8524 


0.9460 


32 


1 


5.0 


1.0138 


0.9055 


0.9715 


84 


3 


6.0 


1.1859 


0.9792 


0.9515 


106 


5 


Z=0.0001 


2.5 


1.1489 


0.7325 


0.8207 


31 


1 


3.0 


1.2758 


0.8159 


0.9691 


40 


1 


3.5 


1.2627 


0.8474 


0.9505 


59 


1 


4.0 


1.1767 


0.8694 


0.9188 


76 


1 


5.0 (VW93) 


1.0374 


0.9340 


0.9297 


138 


2 


5.0 (R75) 


1.2558 


0.9244 


0.9312 


69 


1 



Table 5: Stellar yields for 3M , Z = 0.02 model. 



Isotope i 


A 




yield 




mass(i)i ost 


mass(i)o 






X0(i) 


/ 


Q a 


1 


1 


0600706E- 


■07 


1 


.2625863E-05 


1.2519856E-05 


5 


4468778E- 


-06 


5.4014972E-06 


3.6334887E-03 


p 


1 


-8 


0847383E- 


■02 




1.5124202 


1.5932676 


6 


5246773E- 


-01 


6.8739051E-01 


-2.2644492E-02 


3 He 


3 


2 


2212508E- 


■04 


2 


.2214468E- 


-04 


1.9598588E-08 


9 


5834628E- 

. KJ KJ KJ 7t KJ K\J 1 i 


-05 


8.4555056E-09 


4.0543828 


4 He 


4 


6 


1 267972E- 


■02 


7 


.4012238E- 


-01 


6.7885441E-01 


3 


1 929350E- 

* 1 -7 — -7 *J J V 7 l 1 


-01 


2.9288119E-01 


3.7498605E-02 


7 Li 


7 


-1 


8814418E- 


■08 


6 


.0166285E- 


-09 


2.4831047E-08 


2 


59561 19E- 


-09 


1.0712970E-08 


-6.1567014E-01 


12 C 


12 


1 


4943778E- 


■02 


2 


.2845197E- 


02 


7.9014199E-03 


9 


8555638E- 

- U *J \J *J \J kj KJ J—/ 


-03 


3.4089447E-03 


4.6106154E-01 


13 C 


13 


1 


1 727007E- 


■04 


2 


.1251373E- 


-04 


9.5243653E-05 


9 


1679773E- 


-05 


4.1091393E-05 


3.4852266E-01 


14 N 


14 


3 


3946566E- 


■03 


5 


8388105E- 


-03 


2.4441539E-03 


2 


5189000E- 

. ' J J- kj \J\j\j\j j_j 


-03 


1.0544922E-03 


3.7816754E-01 


15 N 


15 


-4 


7271 697E- 


■06 


4 


.8713855E- 


06 


9.5985552E-06 


2. 


1015467E- 


-06 


4.1411477E-06 


-2.9458165E-01 


16Q 


16 


-1 


2392532E- 


■03 


2 


.1018269E- 


-02 


2.2257522E-02 


9 


06741 53E- 

• \J\J t i i J * } \ j 


-03 


9.6026622E-03 


-2.4908133E-02 


17 


17 


5 


1704166E- 


■05 


6 


0690640E- 


-05 


8.9864743E-06 


2 


.6182328E- 


-05 


3.8770745E-06 


8.2950413E-01 


18Q 


18 


-1 


3407498E- 


■05 


3 


6671925E- 


-05 


5.0079423E-05 


1 


5820502E- 

. 7 ^ 7 — V 7 7 V 7 1 j 


-05 


2.1605989E-05 


-1.3535389E-01 


19p 


19 


3 


871 5425E- 


■06 


4 


.9463952E- 


06 


1.0748529E-06 


2 


1 339065E- 


■06 


4.6372855E-07 


6.6291153E-01 


20 Ne 


20 


-1 


1 501 834E- 


■06 


4 


1906834E- 


-03 


4.1918335E-03 


1 


8078875E- 


-03 


1.8085014E-03 


-1.4747151E-04 


21 Ne 


21 


4 


0954365E- 


■07 


1 


1120977E- 


-05 


1.0711433E-05 


4 


7976600E- 


-06 


4.6212808E-06 


1.6267121E-02 


22 Ne 


22 


1 


990781 4E- 


03 


2 


.3274263E- 


-03 


3.3664503E-04 


1 


0040666E- 


-03 


1.4524026E-04 


8.3967549E-01 


23 Na 


23 


7 


5442040E- 


■05 


1 


.6396206E- 


-04 


8.8520021E-05 


7. 


0734277E- 


-05 


3.8190585E-05 


2.6767361E-01 


24 Mg 


24 


-8 


975621 3E- 


■07 


1 


.3679401E- 


-03 


1.3688377E-03 


5 


9013808E- 

.Jul uuuuu 


-04 


5.9056369E-04 


-3.1310008E-04 


25 Mg 

O 


25 


2 


4309760E- 


■05 


2 


.0399288E- 


-04 


1.7968312E-04 


8. 


8003828E- 


-05 


7.7521494E-05 


5.5079415E-02 


26 Mg 


26 


1 


1070690E- 


■05 


2 


1723099E- 


-04 


2.0616030E-04 


9 


3714836E- 

* KJ 1 1 1 77 ' 7 V 1 ^ 


-05 


8.8944660E-05 


2.2688469E-02 


26 AI 


26 


1 


1710374E- 


■07 


1 


.1710374E- 


-07 


0.0000000E+00 


5 


0519297E- 


-08 


0.0000000E+00 


0.0000000E+00 


27 AI 


27 


1 


2958917E- 


■06 


1 


.5529012E- 


-04 


1.5399423E-04 


6 


6993147E- 

• V -7 -7 '.7 1 1 f J_J 


-05 


6.6438413E-05 


3.6111036E-03 


28 Si 


28 


-1 


.6319100E- 


06 


1 


.7316726E- 


03 


1.7333045E-03 


7. 


.4705458E- 


■04 


7.4780727E-04 


-4.3735650E-04 


29 Si 


29 


1 


.5787009E- 


06 


9 


.2347655E- 


-05 


9.0768954E-05 


3 


.9839368E- 


■05 


3.9160856E-05 


7.4602445E-03 


30 Si 


30 


2 


.4255860E- 


■07 


6 


.2702733E- 


■05 


6.2460174E-05 


2 


.7050359E- 


■05 


2.6947471E-05 


1.6549942E-03 


31p 


31 


2 


0692660E- 

• V M7 ' / w V7 \J WJ-7 


06 


2 


.3703114E- 


■05 


2.1633849E-05 


1 


.0225674E- 


-05 


9.3335875E-06 


3.9643381E-02 


32 s 


32 


-2 


1749875E- 


■06 


1 


.0513477E- 


■03 


1.0535227E-03 


4. 


5355811E- 

* KJ KJ KJ KJ KJ _1_ 1 1 


-04 


4.5452599E-04 


-9.2578813E-04 


33 s 


33 


4 


.0168288E- 


-07 


8 


.9708283E- 


■06 


8.5691454E-06 


3 


.8700723E- 


-06 


3.6970248E-06 


1.9866738E-02 


34g6 


34 


-1 


.7065395E- 


-07 


4 


9333787E- 


-05 


4.9504441E-05 


2 


1282909E- 


-05 


2.1357921E-05 


-1.5279896E-03 


56p e 


56 


-2 


.1018088E- 


-05 


3 


.0820514E- 


-03 


3.1030695E-03 


1 


.3296166E- 


-03 


1.3387711E-03 


-2.9798900E-03 


57 Fe 


57 


1 


.3097335E- 


■05 


8 


.8821595E- 


-05 


7.5724260E-05 


3 


8318201E- 


-05 


3.2670057E-05 


6.79255233E-02 


58p e 


58 


7 


.2211260E- 


-06 


1 


.7037621E- 


-05 


9.8164946E-06 


7. 


.3501383E- 


-06 


4.2351739E-06 


2.3942426E-01 


60p e 


60 


3 


.9827793E- 


-08 


3 


.9827793E- 


-08 


0.0000000E+00 


1 


.7181964E- 


-08 


0.0000000E+00 


0.0000000E+00 


59 Co 


59 


1 


.4924935E- 


-06 


1 


.0411269E- 


-05 


8.9187752E-06 


4 


.4914877E- 


-06 


3.8478670E-06 


6.7170151E-02 


58 Ni 


58 


-2 


.2356980E- 


-06 


1 


.2900950E- 


-04 


1.3124519E-04 


5 


.5655517E- 


-05 


5.6623696E-05 


-7.4899960E-03 


60 Ni 


60 


7 


.0899841E- 


■07 


5 


.2735955E- 


-05 


5.2026957E-05 


2. 


.2750626E- 


-05 


2.2446222E-05 


5.8501000E-03 


61 Ni 


61 


8 


.3594359E- 


-06 


1 


.0636427E- 


-05 


2.2769907E-06 


4 


.5886222E- 


-06 


9.8237228E-07 


6.6940618E-01 


62 Ni 


62 


-7 


.3240267E- 


-06 


5 


.0168190E- 


-08 


7.3741949E-06 


2 


1642876E- 


-08 


3.1814818E-06 


-2.1673145E+00 



•a 
c 



a g represents the abundance of species from 64 Ni to Bi; an increase in g indicates that neutron-captures have occurred beyond the end of the network. 

b 34 S includes the abundance of species between 34 S and Mn. 
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Table 6: PN yields 3M Q , Z = 0.02 model. The first entry is the total mass ejected during the final two 
TPs, 1.502 M in this case. 



Isotope i 


A 




yield 




mass(i) lost 


X0(i) 








N(i)/N(H) 


g 


1 


9. 


.4299139E- 


■08 


8. 


.2076585E- 


-06 


5.4014972E-06 


5. 


.4642769E- 


-06 


8.3678888E-06 


p 


1 


-5 


.1648498E- 


-02 


9 


.8085165E- 


■01 


6.8739051E-01 


6, 


.5300536E- 


■01 


1.0000000E+00 


3 He 


3 


1. 


.4400885E- 


■04 


1 


.4402156E- 


■04 


8.4555056E-09 


9. 


5882846E- 


-05 


4 


.8944388E- 


-05 


4 He 


4 


3 


.9245605E- 


-02 


4 


.7917002E- 


■01 


2.9288119E-01 


3 


1900910E- 


■01 


1 


.2213111E- 


-01 


7 Li 


7 


-1 


.2191149E- 


-08 


3 


.9003480E- 


-09 


1.0712970E-08 


2. 


.5966702E- 


-09 


5 


.6807015E- 


■10 


12 C 


12 


9. 


.4264597E- 


■03 


1 


.4546891E- 


02 


3.4089447E-03 


9. 


.6846428E- 


-03 


1 


.2359064E- 


-03 


13 C 


13 


7. 


.6098098E- 


-05 


1 


.3781973E- 


■04 


4.1091393E-05 


9. 


1753966E- 


-05 


1 


.0808482E- 


-05 


U N 


14 


2 


.2036042E- 


-03 


3 


.7875122E- 


■03 


1.0544922E-03 


2. 


.5215494E- 


-03 


2 


.7581805E- 


■04 




15 


-3 


.0603858E- 


-06 


3 


.1598568E- 


-06 


4.1411477E-06 


2. 


1036856E- 


■06 


2 


.1476960E- 


■07 


16Q 


16 


-7 


.9538208E- 


-04 


1 


.3628369E- 


■02 


9.6026622E-03 


9. 


0731336E- 


-03 


8 


.6840155E- 


■04 


17 


17 


3 


.3530909E- 


■05 


3 


.9354498E- 


■05 


3.8770745E-06 


2. 


6200392E- 


-05 


2 


.3601635E- 


-06 


18Q 


18 


-8 


.6711771E- 


-06 


2 


.3782261E- 


■05 


2.1605989E-05 


1. 


5833122E- 


-05 


1 


.3470301E- 


-06 


19p 


19 


2 


.4399674E- 


-06 


3 


.1365144E- 


-06 


4.6372855E-07 


2. 


.0881453E- 


■06 


1 


.6830241E- 


■07 


20 Ne 


20 


-8 


.7078661E- 


-07 


2 


.7156025E- 


■03 


1.8085014E-03 


1. 


.8079218E- 


-03 


1 


.3843084E- 


■04 


21 Ne 


21 


2 


.5139525E- 


■07 


7. 


.1928253E- 


-06 


4.6212808E-06 


4. 


7886483E- 


■06 


3 


.4920211E- 


■07 


22 Ne 


22 


1. 


.2537986E- 


-03 


1 


.4719579E- 


■03 


1.4524026E-04 


9. 


7996101E- 


-04 


6 


.8213347E- 


•05 


23 Na 


23 


4 


.7911446E- 


-05 


1 


.0527590E- 


■04 


3.8190585E-05 


7. 


0087794E- 


-05 


4 


.6665705E- 


-06 


24 Mg 


24 


-5 


.7858415E- 


-07 


8. 


.8648207E- 


■04 


5.9056369E-04 


5, 


.9017848E- 


-04 


3 


.7657839E- 


-05 


25 Mg 


25 


1 


.4788740E- 


■05 


1 


.3123048E- 


■04 


7.7521494E-05 


8. 


.7367145E- 


-05 


5 


.3516951E- 


-06 


26 Mg 


26 


6 


.7903602E- 


-06 


1 


.4039036E- 


■04 


8.8944660E-05 


9. 


.3465365E- 


-05 


5 


.5050414E- 


-06 


26 A1 


26 


7 


.8603932E- 


-08 


7. 


.8603932E- 


-08 


0.0000000E+00 


5. 


.2330840E- 


-08 


3 


0822480E- 


-09 


27 Al 


27 


8 


.2996121E- 


■07 


1 


.0062428E- 


■04 


6.6438413E-05 


6. 


6990964E- 


-05 


3 


.7995808E- 


-06 


28 Si 


28 


-1 


0048971 E- 


-06 


1 


1 222448E- 

. J_ Lj Lj ilj t: t: O 1 I 


-03 


7.4780727E-04 


7. 


471 3834E- 


-04 


4 


.0862618E- 


-05 


29 Si 


29 


9 


.7531665E- 


■07 


5 


.9797174E- 


■05 


3.9160856E-05 


3 


.9810177E- 


-05 


2 


.1022256E- 


-06 


30 Si 


30 


1 


.4533725E- 


■07 


4 


.0621992E- 


■05 


2.6947471E-05 


2. 


.7044231E- 


-05 


1 


.3805007E- 


-06 


31p 


31 


1 


.2756063E- 


-06 


1 


.5295192E- 


■05 


9.3335875E-06 


1. 


.0182827E- 


-05 


5 


.0302543E- 


■07 


32 g 


32 


-1 


.3433164E- 


-06 


6 


.8138086E- 


■04 


4.5452599E-04 


4. 


5363166E- 


-04 


2 


.1708838E- 


-05 


33 g 


33 


2 


5333611E- 


■07 


5 


.8064797E- 


■06 


3.6970248E-06 


3 


8656840E- 


■06 


1 


.7938892E- 


■07 


34 g 


34 


-1 


.0619988E- 


-07 


3 


.1974625E- 


■05 


2.1357921E-05 


2. 


.1287216E- 


-05 


9 


.5878931E- 


■07 


56p e 


56 


-1 


.3052253E- 


-05 


1 


.9978590E- 


■03 


1.3387711E-03 


1. 


.3300815E- 


■03 


3 


.6372523E- 


-05 


57 Fe 


57 


8 


.2095721E- 


-06 


5 


.7281879E- 


■05 


3.2670057E-05 


3 


.8135608E- 


-05 


1 


.0245640E- 


-06 


58p e 


58 


4 


.4331669E- 


-06 


1 


.0794642E- 


■05 


4.2351739E-06 


7. 


1865697E- 


■06 


1 


.8974787E- 


■07 


60p e 


60 


2 


.0917762E- 


-08 


2 


.0917762E- 


-08 


0.0000000E+00 


1. 


.3926072E- 


-08 


3 


.5543535E- 


■10 


59 Co 


59 


9 


.3833069E- 


■07 


6 


.7180481E- 


-06 


3.8478670E-06 


4, 


.4725639E- 


■06 


1 


.1608812E- 


■07 


58 Ni 


58 


-1 


.3969257E- 


-06 


8. 


.3655119E- 


■05 


5.6623696E-05 


5. 


.5693683E- 


-05 


1 


.4704869E- 


-06 


60 Ni 


60 


4 


.3515684E- 


■07 


3 


.4150675E- 


■05 


2.2446222E-05 


2 


.2735929E- 


-05 


5 


.8028951E- 


■07 


61 Ni 


61 


5 


.3855811E- 


■06 


6 


.8611607E- 


■06 


9.8237228E-07 


4. 


.5678416E- 


■06 


1 


1467385E- 


■07 


62 Ni 


62 


-4 


. 75435 13E- 


-06 


2 


.4417719E- 


■08 


3.1814818E-06 


1 


.6256180E- 


-08 


4 


.0152268E- 


-10 
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